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We investigate the presence of off-diagonal long-range order in a harmonically confined two- 
' dimensional Bose gas. In the noninteracting case, an analytical calculation of the the finite- 

temperature one-particle density martix provides an exact description of the spatial correlations 
known to be associated with the existence of a Bose-Einstein condensate below the transition tem- 
£N) ' perature Tj ' . We treat the effects of repulsive interactions within the semiclassical Hartree-Fock- 

, | | Bogliubov approximation and find that even though the system remains in the same uncondensed 

O , phase for all T > 0, there appears to be a revival of off-diagonal long-range order for temperatures 

is not the BEC state. 
PACS numbers: 03.75.Fi,05.30. Jp 
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T < tJ ' . We suggest that this reentrant order is related to a phase transition in the system which 



I. INTRODUCTION 



q_) ' Recent advances in the controlled fabrication and manipulation of ultra-low temperature trapped Bose gases |l], |^| 
has finally made it possible to experimentally investigate the existence of off-diagonal long-range order (ODLRO), 
i.e., Bose-Einstein condensation (BEC), in dimensions lower than three. This remarkable experimental achievment 
has rekindled an interest in the classic theoretical problem concerning the existence of BEC in two-dimensions (2D). 
Although it is well known that finite-temperature BEC can never occur in a homogeneous 2D Bose gas || , the ideal 

trapped Bose gas can undergo BEC below some critical temperature Tc°\ In the case of a 2D isotropic harmonic 
trap with confining frequency ojq (unless stated otherwise, we always implicitly assume this confinement geometry), 

Bagnato and Kleppner have shown that \/pc = huJo^/GN/?: 2 [Q, where (3 — l/(kgT); Shevchenko later extened this 
result to include more general power law potentials [H. In the thermodnyamic limit, viz., N — > oo, ujq — > 0, such that 
N^-^luq = constant, BEC is also known to theoretically occur with the critical temperature remaining unchanged. 
' Note that this is not the usual thermodynamic limit, which in 2D would demand that Nujq — constant, resulting in 
no BEC (see also Ref. [|| for an alternative thermodynamic limit in the case of an interacting Bose gas). 

If repulsive interactions between bosons are included, (we only consider repulsive interactions in this paper), then a 
\& [ definitive answer to the question of BEC in the trapped 2D gas is not so clear. In fact, there is no a priori reason for 
the existence of BEC in a confined two or three-dimensional system — BEC is a purely kinematical phenomenon in the 
sense that even the ideal Bose gas can undergo the phase transition. In contrast, a transition to the superfluid state 
is wholey dynamical and cannot occur without interactions. Therefore, one could sensibly assume that interactions 
prevent the spatial redistribution of atoms required to achieve Bose gas degeneracy, and therefore destroy the BEC 
phase transition altogether. 

It is only recently that a rigorous proof has been given to show that there is no finite temperature BEC in the 
interacting trapped 2D Bose gas in the thermodynamic limit 0. Of course, in actual experiments, the theormodyamic 
I ■ limit is never fully realized, and it becomes necessary to consider systems consisting of a finite number of atoms confined 
in a trap. In this regard, Lieb and Seiringer || have proved that the finite, interacting trapped Bose gas does have 
100% BEC at T = 0. To date however, no irrefutable case has been made (experimentally or theoretically) in support 
of finite temperature BEC in an interacting finite 2D system. 

In a recent paper, Petrov et al. || have suggested that that there are in fact two BEC regimes for finite, trapped, 
(quasi)-2D systems. The general argument supporting their claim hinges on an approximate analytical calculation 
of the one-particle density martix from which information about the phase correlations of the condensate can be 
obtained. For large particle number and sufficiently low temperatures, Petrov et al. deduce that the phase fluctuations 
are governed by ~ Tln(iV) Jl. Thus, at temperatures well below T c , fluctuations in the phase are supressed and one 
has a "true condensate" whereas for T < T c , the phase fluctuations are enhanced and one has a "quasi-condensate" . 
Similar results have recentl y b een obtained by Bogliubov et al., in strictly two-dimenions B. 

However, Bhaduri et al. |lfl have suggested that there is no BEC phase transition in a finite, interacting 2D Bose 
gas (see also Refs. O, O, [13| for similar conclusions). In their investigation, a semiclassical approximation to the 
Hartree-Fock-Bogliubov (HFB) mean-field theory was used to investigate the thermodynamics of the trapped 2D Bose 
gas. Although the semiclassical HFB approach could not shed any light on the question of phase fluctuations, it was 
shown that the HFB equations could be solved self-consistently all the way down to T = provided one does not 
invoke the presence of a condensate. This finding seems to indicate that that if no a priori assumption of a condensate 
is made (as it should be), the system remains in the same uncondensed phase at all temperatures. (Note that applying 
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the same HFB approach in 3D requires a condensate in order to obtain self-consistent solutions below Tc°\) In spite 
of the absence of a condensate order-parameter however, the density distribution for the uncondensed phase was found 

to clearly exhibit a large spatial accumulation of atoms near the center of the trap for T < T c , similar to what is 
found in systems with a condensate. This suggests that there may be some sort of "revival" of ODLRO in the system, 
even though it does not lead to the formation of a condensate. Since the off-diagonal one-particle density matrix can 
provide information about phase correlations and ODLRO in the system, an investigation of the off-diagonal density 
martix within the HFB approximation needs to be performed. This is indeed the central motivation for the present 
work. 

The rest of our paper is organized as follows. In Sec. II we investigte the noninteracting 2D gas in a harmonic trap 



where an exact expression for the finite temperature correlation function can be obtained. In Sec. HI, we calculate 
the correlation function for the interacting system within the semiclassical HFB approximation and investigate its 
temperature dependence both above and below Tc . Finally, in Sec. ^v| , we present our concluding remarks. 

II. NONINTERACTING GAS 

In this section, we wish to pay special attention to the off-diagonal, one-particle density matrix (r, r'; (3) of the 
noninteracting, confined 2D Bose gas. Our motivation lies in the fact that ODLRO in the one-particle density matrix 
and BEC are interrelated Jl4| (see below for details). For the isotropic HO, p^(r, r'; j3) can be calculated analytically, 
and we have an exact result from which we can quantitatively test the validity of the semiclassical approximation. 
We will be particularly interested how well the seimclassical approximation describes the phase correlations of the gas 
below the BEC critical temperature Tc°\ To the best of our knowledge, this is the first such study in two-dimensions. 

By definition, the off-diagonal one-particle density matrix is given by |l5| ] 



p(°)(r,r';/3)=^^f(r)^(r') t nf , (1) 



where V>j (r) are exact eigenfunctions of the <i-dimensional harmonic oscillator (HO) and rSp is the Bose distribution 
function. It is most convenient to evaluate ([j]) for the ID case, and then use the separability of the HO potential to 
obtain the general result. Explicitly, we have 



- ' , (2 ) 



n=0 



where e„ are the energy levels of the HO, p is the chemical potential, and we have scaled all energies and lengths 
by hojQ and £g — y^h/muio, respectively. We now set the zero of energy to correspond to the lowest trap level, viz., 
ri = 0, and obtain from (||) 

P W(x,x';P) = -L f; J2 e-^ 2+ *' 2 > Hn{X ^ {X ' ] e-^ n , (3) 



3=1 n=0 



where we have used [hi 



E 



1 



e 



/3(e:„-A0 - 1 



(4) 



The summation over trap levels can be perfomed exactly Jl^j , and we finally obtain in any dimension 



p {0) (r,r';P) = ^ 9 C~ 1/4[|r+r ' 1 tanh(. 7) 3/2) + |i-r'| 2 coth(j/3/2)] _ ^ 

.7=1 t 71 "^ 



Putting d = 3 in the above equation reproduces the result obtained by Barnett et al, The chemical potential is 
obtained by enforcing that the total particle number N, remain fixed, viz., 
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FIG. 1: Exact quantum-mechanical densities for temperatures (starting from the lowest curve) T/Tr = 1.2, 1.0,0.9,0.8,0.6. 
The figure inset shows the fractional occupancy No/N of the lowest oscillator level for N — 10 4 atoms (solid curve) and the 
thermodyanmic limit (dashed line) as given by Eq. (7). 



The diagonal part of the one-particle density matrix yields the exact quantum-mechanical finite temperature density, 
p(°'(r), of the trapped gas. The densities for various temperatures are shown in Fig. 1 for d = 2 and N = 10 4 atoms. 
Figure 1 illustrates that even though we are not in the thermodynamic limit, the single-particle density distribution 
exhibits an enormous accumulation of atoms near the center of the trap for T < Tj . This type of behaviour is 
consistent with the confined gas having undergone a phase transition to the BEC state. The usual approach for 
confirming that this phase transition is indeed BEC is to examine the temperature dependence of the fractional 
occupancy of the lowest oscillator level, No/N. The inset to Fig. 1 shows that below T/Tc ~ 1, the ground state 
of the system becomes macroscopically occupied, with already 50% of the atoms in the lowest level at T/Tc ~ 0.7. 
Remarkably, even at iV = 10 4 atoms, the fractional occupancy closely follows the quadratic temperature dependence 
obtained strictly in the thermodynamic limit (dashed line in figure inset), viz., 

^ = l-(— V (7) 

Although N /N gives us a strong indication that the finite system has a BEC transition at T/T^ ~ 1, it docs not 
provide us with a definitive test for the presence of ODLRO in the system. 

It is well known since the early work of Penrose and Onsager | fl4| that for a homogeneous system, the phenomenon 
of BEC is intimately related to the presence of ODLRO in the off-diagonal one-particle density matrix, namely, 

lim p(°)(r,r';/3) ^ . (8) 
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In fact, the low temperature behaviour of p^(r, r') is directly related to the phase fluctuations of the condensate O] 
and the non-zero limiting value of Eq. (8) yields the value of the order parameter (i.e., condensate fraction) in the 
system. For finite inhomogeneous systems however, the above criterion has to be modified to account for the fact 
that ODLRO cannot extend beyond the physical size of the system. This is accomplished by defining a normalized 
correlation function (sometimes called the first-order coherence function) by 

7 (°W-m = : ? (0) (r,r';/3) (9) 

V ' VP (0) (r,r;/3)Vp (0) (r',r';/3) ' 

such that 

lim 5 <°>(r,r';/3)^0. (10) 

The above correlation function describes a local measure of coherence in the system. For example, the ability of two 
BEC's, intially separated by |r — r'|, to form an interference pattern is quantified by g^(r, r') |jl9). In this sense, the 
correlation function also characterizes local fluctuations of the phase of the condensate order-parameter. The quantity 
£ is the large-distance length-scale over which the ODLRO of the finite system is manifested. At zero temperature, 
£ = £q, whereas for T ^ 0, £ is taken to be the spatial extent of the single-particle density p^^r). To simplify 
our calculations, we fix one of the coordinate, r' to be at the center of the trap. Note that in Eq. (9), dividing the 
off-diagonal density matrix by the square root of the single-particle densities removes any effects of the local-density 
of particles on the correlation function, and gives g^°> (r, r; (3) — 1 in the case of perfect coherence. 

The solid curves in Fig. 2 depict the correlation function g^°\r, 0) for d = 2 and N = 10 4 at various temperatures. 
From left to right, we have T/Tc — 1.0,0.9,0.8,0.6,0.4.0.2,0.01. Above the critical temperature, the atoms remain 
in the uncondensed phase and the system shows no ODLRO. However, as the gas is cooled below the transition tem- 
perature, we see an abrupt jump in the large-distance decay of the spatial correlations. A direct application of Eq. (10) 

yields condensate fractions of N /N = 0.02,0.47,0.71,0.84,0.93, for temperatures T/T c (0) = 1.0,0.8,0.6,0.4,0.2, re- 
spectively. Thus, below Tc°\ the finite noninteracting 2D Bose gas contains a non-zero order-parameter (i.e., a 
BEC). 

It also proves very useful to compare the spatial extent of the single-particle density, p^ (r) , to the spatial decay of 
the correlation function </°)(r,0). This information is also presented in Fig. 2 where the exact quantum mechanical 
densities are represented by dashed lines. These densities have been scaled by their value at r = to allow for a 
comparison of the length scales between p(°)(r) and g(°\r,0). From right to left, the dashed lines correspond to 

densities at T/Tc — 1.0,0.9,0.8,0.6. Densities at lower temperatures are almost indistinguishable on the scale of 
the figure. We note that at low temperatures, the correlation function is practically a constant g^°'(r, 0) = 1 over 
the entire region for which the density (i.e., physical size of the system) of the gas is nonzero. This means that 
phase fluctuations are very small, and the density distribution corresponds to that of a "true condensate" . At higher 
temperatures (i.e., 0.9 < T/T^ < 1.0) the single-paricle density can become comparable to that of </°)(r, 0), and the 
gas can be said to be a "quasi-condensate" in the sense that the phase correlations do not extend to all of the atoms 
in the trap. For T/Tc > 1-0, the spatial extent of the single-particle density greatly exceeds characteristic length of 
the correlations, and system is in the normal state. 

Having obtained the exact expression for p( )(r,r';/3), we now wish to examine its semiclassical approximation, 
namely, its behaviour when fc^T ^ Tiujq. It is easy to show that Eq (5) reduces to 

p(0)(ry) = (^^E^^pH^^^ (ii) 



1 



5 rf/2 (exp [f3(p - [V(r) + V{v')]/2)] , exp [-|r - r'| 2 )/2/3] ) , (12) 



(27r/3) rf / 2; 

where V(r) = l/2r 2 , and as in Ref. P0| we have introduced the generalized Bose function 



oo j l/j 

9a(x,y) = ^2^j^- (13) 



For r = r' and d = 2, Eq. (12) becomes 



p(°)(r) 



1 [ d 2 p 



(2nf 



exp 



(£ + ir* + - M ) P 



(14) 
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FIG. 2: Exact off-diagonal correlation function (solid curves) for temperatures (from left to right) T/T^ = 
1.0,0.9,0.8,0.4,0.2,0.01. The dashed lines correspond to the scaled exact quantum mechanical densities at temperatures 
(from right to left) T/Tc = 1.0,0.9,0.8,0.6 The filled circles are the semiclassical approximation to the off-diagonal density 
matrix as given by Eq. (16). 



with the normalization condition 



N 



(15) 



If one attempts to solve Eqs. (14,15) for all T, it is readily shown that there exists a temperature T* below which the 
equations can no longer be satisfied p0[ O. In fact, the temperature T* is identical to the the critical temperature 
obtained by Bagnato and Kleppner ||. Thus, in order to obtain solutions below Tc°\ a condensate order parameter is 
required. We introduce the macroscopic occupation of the ground state in the one-particle density matrix according 
to the ansatz 



/)(r,r') = JVoC J (r)C , (r') + 



(0)/ 



1 



rff d/ a(exp[/30*-[V(r)+V(r / )]/2)],exp[-|r-r / | 2 /2/3]) , (16) 



where ^q "* (r) is the zero temperature ground-state eigenfunction of the HO and 



Nn = 



The 2D density is then given by Eq. (16) with r = r', 



p(°>(r)=iV |CW + 



(0), 



(2tt)2 



d 2 p 



exp 



- 1 



(17) 



(18) 
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The semiclassical densities obtained from Eq. (18) are indistinguishable from the exact densities obtained from 
Eq. (5). The semiclassical correlation function, Eq. (16), at various temperatures is represented in Fig. 2 as filled 
circles. The superb agreement between the exact (solid curves) and semiclassical results shows that the semiclassical 
approximation is valid well below the critical temperature provided one takes into account the macroscopic occupation 
of the ground state separately. If the atoms had remained in the uncondensed phase, Eqs. (14,15) would have been 
sufficient to describe the gas for all T > 0. 



III. INTERACTING GAS 



The thermodynamic properties of the interacting Bose gas are most commonly described by the self-consistent 
mean- field HFB equations. In most calculations (see Ref. [gjj ^2) and references therein for details), the semiclassical 
approximation is applied to the discrete set of coupled Bogliubov equations, leading to the considerably simplified 
semiclassical HFB theory. The semiclassical HFB approximation is well-known to be an excellent tool for describing 
the interacting 3D Bose gas [^l], [22] ]. The 2D version of the semiclassical HFB theory is formally identical to the 
3D version and has already been examined by Mullin (^3). We present only the essentials of the model here for 
completeness. 

In the semiclassical HFB approach, the condensed and uncondensed atoms are described by a set of coupled, self- 
consistent equations. The macroscopic condensate wavefunction, ^o(r), is characterized by the finite-temperature 
Gross-Pitaevskii equation which is given by 



1 r 



1 2 

r 



^j + 2 7 p(r) -7Po(r) 



Mr) = 



(19) 



where 7 is the coupling constant associated with the assumed repulsive zero-range two-body pseudo-potential. The 
total density, p(r), is made up of a condensate po( r ) = l - 0o( r )| 2 I and a noncondensate density, 



p T (r) = J { [u 2 (p, r) + v 2 {p, r)] /(p, r) + v 2 {p, r)} d 2 p 



where 



(20) 



/(P,r) 



(21) 



and 



The total number of particles satisfies 



u 2 (r,p) = 

v 2 {r,p) = 
p 2 1 



A + e 

A-E 



2e 



A = y + -r 2 - fi + 2 7 p(r) 



:(p,r) = VA 2 - ( 7( o (r)) s 



N = J [p (r) + pt(t)] d 2 r. 
In the absence of a condensate, it is easy to show that the above set of HFB equations reduce to 

1 r d 2 P 



Pr(r) 



2n 



exp 



(£ + |r2 + 2 7PT (r) - M ) (3 



1 



(22) 
(23) 
(24) 

(25) 

(26) 

(27) 
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FIG. 3: Semiclassical densities for N — 10 4 atoms with 7 = 0.16 (solid curves) and 7 = 0.31 (dashed curves). The lowest lying 
curves correspond to T/tJ ' = 1.1 and the higher lying curves to T/Tj ' = 0.8. The figure inset shows the chemical potential 
as a function of temperature. The long-dashed-short-dashed curve corresponds to 7 = 0. The value 7 = 0.31 is consistent with 
an effective 2D interaction strength for Rb 87 Jicj]. 



with 

N = J p T {r) d 2 r. (28) 

Equations (27,28) are equivalent to the so-called self-consistent Thomas- Fermi approximation used in Refs. jL0[ [Til . 
Putting 7 — in Eq. (27) yields an expression that is identical to the semiclassical approximation of the exact diagonal 

density matrix, viz., Eq. (14), found in the previous section; recall from Sec. II that below Tc°\ Eqs. (27,29) cannot 
be solved self-consistently. However, if one puts any finite 7 > in Eq. (27), self-consistent solutions to the HFB 
equations can be obtained all the way down to T = [|l0| [0|, [l^]. To illustrate this, we show in Fig. 3 the density and 
chemical potential (figure inset) of the gas for 7 = (long-dash-short-dashcd line), 7 = 0.16 (solid line), and 7 = 0.31 
(dashed line). In contrast to the noninteracting 2D gas (where a condensate order-parameter has to be introduced 

below Tc )j the interacting system remains in the uncondensed phase for all T > 0. 

Although no condensation occurs, it is interesting to compare the densities of the noninteracting (Fig. 1) and 

interacting gases (Fig. 3) just above and below Tc°\ The essential point to be made from this comparison is that 
in both cases, there is an increase in the density of atoms near the center of the trap for T < Tc°\ While the 
noninteracting density shows a much larger enhancement at r = 0, the interacting gas also displays a similar (albeit 
diminished) behaviour. Since an increase in the density of atoms near the center of the trap is associated with 
ODLRO in the noninteracting gas, it is worthwhile investigating the off-diagonal density matrix for the interacting 
2D gas within the semiclassical HFB approximation. 
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The 2D semiclassical HFB off-diagonal density matrix can be obtained from the leading order term in the Wigner- 
Kirkwood semiclassical expansion of the quantum mechanical off-diagnonal density matrix (24). The noncondensed 
particles are treated as bosons in an effective potential, and we obtain 

P(r, r'; 0) = ^5i(exp - [V cS (r) + V eS (r')]/2)] , exp [-|r - r'| 2 /2/3] ) , (29) 

where 

Veu(r) = ir 2 + 2 7 p(r) . (30) 

Since there is no condensate for 7 > 0, we have not included the first term in Eq. (16). Of course, for r = r', Eq. (29) 
reduces to Eq. (27). 

As in Sec. II, the normalized correlation function for the finite sized system is given by 



where again, we have fixed r' = 0. The temperature dependence of g{v, 0; (3) is shown in Figs. 4, 5 for coupling 

strengths 7 = 0.16,0.31, respectivley. Above Tc°\ the interacting and noninteracting (see Fig. 2) systems have 

similar correlations. However, as we cool the gas to temperatures just below Tc°\ there is an drastic change in the 
behaviour of the interacting off-diagonal density matrix, which to the best of our knowledge, has never been noticed 
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FIG. 5: As in Fig. 4 but for 7 = 0.31. 
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before. In particluar, we observe that the decay of the correlations in the interacting system are qualitatively different 
from the noninteracting gas. For the interacting gas, values of r > Iq lead to the tail of the correlation function 
decreasing slowler than the noninteracting case. In contrast, for r < £0, the interacting correlation function has a 
much faster spatial decay than the noninteracting system. We have already discussed the fact that the near constancy 
of the noninteracting correlation function, g(°\r,r'), over the size of the sample (r ~ £0) is indicative of negligible 

phase fluctuations in the condensate Jm). Thus, the re-appearance of ODLRO in g(r, r') below suggests that there 
may be larger fluctuations in the phase as compared to the noninteracting gas. If we apply Eq. (10) to g(r, r') (with 
7 = 0.16), we obtain 0.002,0.25,0.58,0.70 at T/T c (0) = 1.0,0.8,0.6,0.4 respectively, for the asymptotic behaviour of 
the correlation fucntion. Even though we may not be able to interpret these numbers as condensate fractions, it is 
tempting to conjecture that they represent a nonzero value of an order-parameter in the system. 

We can be more quantitative about the nature of the correlations displayed by g(r, r') below Tj ^ as follows. In the 
interacting homogeneous 2D Bose gas, Popov [B5| has already calculated the low temperature correlation function: 



/(r,r';/3)«i?r« , 



where 



R : 



2tt/3 Ps 



(32) 
(33) 



and p s is the so-called superfiuid density. The algebraic ODLRO exhibited by g(r, r') implies that at low temperature, 
correlations can extend over macroscopic (i.e., size of the system) distances, although phase fluctutations will be 
nonzero. This is reminiscent of the quasi-condensate described in Ref. [SL and it is only identically at T = that the 
system will attain a "true condensate" . 
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If the system is finite (i.e., in a 2DHO potential), then Bogliubov et al. 0, have shown that 

_ |r-r'| 1 

R= j^wy a ^^Msj 7 (34) 

where S = (r + r')/2. Eq. (34) is formally identical to (33) except for the slowly varying factor p(S) which is 
the inhomogeneous boson density. This additional spatial dependence results in a deviation from the strict power 
law dependence of the uniform gas. We have found that at low temperatures, the power law dependence in Eq. (34) 
agrees well with the observed numerical behaviour of g(r, r'; (3). This suggests that at low temperatures, the interacting 
system exhibits algebraic ODLRO similar to that of the superfluid state in the uniform Bose gas. At exactly T = 0, 
the gas will have a true condensate. 



IV. CONCLUSIONS 



We have investigated the finite-temperature correlations in both the ideal and interacting trapped, 2D Bose gas. 
Our study of the noninteracting gas provided an exact framework from which we could compare the validity of the 
semiclassical approximation over a wide temperature range. Our results indicate that the semiclassical approach 
works extremetly well, even below the critical temperature, provided the macroscopic occupations of the ground state 
is treated separately. A careful examination of the correlation function reveals that, as in Ref. there are two BEC 
regimes for the noninteracting gas. Namely, for temperatures 0.9 < T/Tj ' < 1.0, the system is a quasi-condensate in 
the sense that the phase correlations are comparable to the spatial extent of the single-particle density. On the other 
hand, when the temperature is lowered to T/Tj ^ < 0.9, phase fluctuations become negligible, and the gas is a true 
condensate. It is straightforward to apply our results to study the finite N dependence on the coherence properties 
of trapped gases in dimensions d = 1, 2, 3. (see e.g., Ref. [Q for work relating to 3D). 

In the interacting gas, we found that even though the system does not contain a condensate, there is a revival 
of ODLRO for T/Tj ^ % 1. The low temperature behaviour of the interacting correlation function was found to 
be consistent with an algebraic (i.e., power law) decay in the phase correlations, similar to what is found in the 
uniform 2D superfluid Bose gas. Whether or not this re-entrant ODLRO is to the superfluid state is presently being 
investigated. 

Finally, it is interesting to note that Margo and Ceperely |^6| have found conclusions analgous to ours for the 
uniform 2D charged Bose fluid (CBF) with ln(r) interactions. Specifically, they found that even though the 2D CBF 
does not undergo a condensation, the noncondensed fluid exhibits a power-law decay of the one-particle density matrix 
similar to what is found in the present work. It is not immediately clear what relevance their findings may have on 
our results for the interacting, inhomogeneous 2D Bose gas. 
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